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O ABSTRACT 

(N ' 

Context. The angular diameter distances toward galaxy clusters can be determined with measurements of Sunyaev- 
Zel'dovich effect and X-ray surface brightness combined with the validity of the distance-duality relation, Dl{z)(1 + 
z) I / Da(z) = 1, where Dl(z) and Da(z) are, respectively, the luminosity and angular diameter distances. This combi- 
nation enables us to probe galaxy cluster physics or even to test the validity of the distance-duality relation itself. 
Aims. We explore these possibilities based on two different, but complementary approaches. Firstly, in order to constrain 
the possible galaxy cluster morphologies, the validity of the distance-duality relation (DD relation) is assumed in the 
ACDM framework (WMAP7). Secondly, by adopting a cosmological-model-independent test, we directly confront the 
angular diameters from galaxy clusters with two supernovae la (SNe la) subsamples (carefully chosen to coincide with 
the cluster positions). The influence of the different SNe la light-curve fitters in the previous analysis are also discussed. 
Methods. We assumed that n is a function of the redshift parametrized by two different relations: n(z) = 1 + rjoz, and 
rj(z) = 1 + r)oz/(l + z), where r/o is a constant parameter quantifying the possible departure from the strict validity 
of the DD relation. In order to determine the probability density function (PDF) of 770, we considered the angular 
Oh! diameter distances from galaxy clusters recently studied by two different groups by assuming elliptical and spherical 

isothermal /3 models and spherical non-isothermal P model. The strict validity of the DD relation will occur only if the 
5^ ' maximum value of 770 PDF is centered on 7/0 = 0. 

Results. For both approaches we find that the elliptical /? model agrees with the distance-duality relation, whereas the 
non-isothermal spherical description is, in the best scenario, only marginally compatible. We find that the two-light 
curve fitters (SALT2 and MLCS2K2) present a statistically significant conflict, and a joint analysis involving the dif- 
ferent approaches suggests that clusters are endowed with an elliptical geometry as previously assumed 
Conclusions. The statistical analysis presented here provides new evidence that the true geometry of clusters is elliptical. 
^ . In principle, it is remarkable that a local property such as the geometry of galaxy clusters might be constrained by a 

CO ' global argument like the one provided by the cosmological distance-duality relation. 
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\ 1. Introduction null geodesies in a Riemannian spacetime (see Ellis 1971, 

• 2007, and references therein). 

The reciprocity law or reciprocity theorem, proved long ago Etherington's reciprocity law can be presented in vari- 

by Etherington (1933), is a fundamental keystone for the ous ways, either in terms of solid angles or relating various 

interpretation of astronomical observations in cosmology, cosmological distances. Its most useful version in the con- 

■ It states that if source and observer are in relative mo- text of astronomical observations, sometimes referred to as 

•jj] I tion, solid angles subtended between the source and the the distance- duality (DD) relation, relates the luminosity 

03 ■ observer are related by geometrical invariants where the distance, D L , with the angular diameter distance, D A , by 

redshift z measured for the source by the observer enters means of the following equation 
in the relation. The core idea of this law comes from the 

invariance of various geometrical properties when there is Dl 

a transposition between the roles of source and observer JJ A 
in astronomical observations. Proofs were presented in the 

context of relativistic geometrical optics, where it comes as Since this result is easil y proved in Friedmann-Lemaitre- 

a consequence of the geodesic deviation equation, as well Robertson- Walker (FLRW) cosmologies, perhaps this is 

as in the context of relativistic kinetic theory, where it is the reason why the generality of the relation above is not 

based on the Liouville integral for collision-free photons. full y appreciated by most authors. Indeed, this law is com- 

The fundamental hypothesis behind the reciprocity law is pletely general, valid for all cosmological models based on 

the one made in General Relativity that light travels along Riemannian geometry, being dependent neither on Einstein 

held equations nor on the nature of matter. Therefore, the 

DD relation is valid for both homogeneous and inhomoge- 

* e-mail: holanda@astro.iag.usp.br neous cosmological models, requiring only that source and 

** e-mail: limajas@astro.iag.usp.br observer are connected by null geodesies in a Riemannian 

*** e-mail: mbr@if.ufrj .br spacetime and that the number of photons is conserved. 
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The DD relation plays an essential role in mod- 
ern cosmology, ranging from gravitational lensing studies 
(Schneidder, Ehlers & Falco 1992) to analyses of galaxy dis- 
tribution and galaxy clusters observations (Lima, Cunha & 
Alcaniz 2003; Cunha, Marassi & Lima 2007; Rangel Lemos 
& Ribeiro 2008; Ribeiro 1992, 1993, 2005; Ribeiro & Stoeger 
2003; Albani et al. 2007; Mantz et al. 2010; Komatsu et 
al. 2011), as well as the plethora of cosmic consequences 
from primary and secondary temperature anisotropies of 
the cosmic microwave blackbody radiation (CMBR) ob- 
servations (Komatsu et al. 2011). Other consequences of 
Etherington's reciprocity relation are the temperature shift 
equation T D = T e /(l+z), where T Q is the observed tempera- 
ture and T e is the emitted temperature, a key result for an- 
alyzing CMBR observations, as well as the optical theorem 
that surface brightness of an extended source does not de- 
pend on the angular diameter distance of the observer from 
the source, an important result for understanding lensing 
brightness (Ellis 2007). In this connection, we also observe 
that any source of attenuation, such as "gray" intergalactic 
dust or exotic photon interaction, contributes to violate the 
DD relation because its proof is based on the conservation 
of the average number of photons. 

Although taken for granted in virtually all analyses in 
cosmology, Eq. ((T|) is in principle testable by means of as- 
tronomical observations (Uzan, Aghanim & Mellier 2004; 
Basset & Kunz 2004; Holanda, Lima & Ribeiro 2010; Li, 
Wu & Yu 2011, Nair, Jhingham & Jain 2011). If one is 
able to find cosmological sources whose intrinsic luminosi- 
ties (standard candles) as well as their intrinsic sizes (stan- 
dard rulers) are known, one can determine both Dl and Da 
and, after measuring the redshifts, test the cosmic version 
of Etherington's result as given by the equation above. Note 
that ideally both quantities must be measured in a way that 
does not utilize any relationship coming from a cosmologi- 
cal model, that is, they must be determined by means of in- 
trinsic astrophysically measured quantities. Therefore, the 
ideal way of observationally testing the reciprocity equa- 
tion ([lj would require independent measurements of intrin- 
sic luminosities and sizes of sources, that is, independent 
determinations of Da and Dl for a given set of sources. 

There are, however, less-than-ideal methods to test Eq. 
(JTJl. These usually assume a cosmological model suggested 
by a set of observations, apply this model in the context 
of some astrophysical effect, thereby trying to see if the 
DD relation remains valid. In this context, Uzan, Aghanim 
& Mellier (2004) argued that the Sunyaev-Zel'dovich effect 
plus X-ray techniques for measuring Da (z) from galaxy 
clusters (Sunyaev & Zel'dovich 1972; Cavaliere & Fusco- 
Fermiano 1978) is strongly dependent on the validity of 
this relation (see details in the next section). Briefly, in 
the context of this phenomenon one may consider the dif- 
ferent electronic density dependencies combined with some 
assumptions about the galaxy cluster morphology in order 
to evaluate its angular diameter distance with basis on Eq. 
((U, such that 
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where Sxo is the central X-ray surface brightness, T e o is 
the central temperature of the intra-cluster medium, A e ijo 
is the central X-ray cooling function of the intra-cluster 
medium, A To is the central decrement temperature, and 8 C 



refers to a characteristic scale of the cluster along the line of 
sight (l.o.s.), whose exact meaning depends on the assump- 
tions adopted to describe the galaxy cluster morphology. 

On the other hand, in order to test the validity of DD 
relation it is convenient to assume a deformed expression: 



and from Eq. (2) we have (see section 4 for details) 



Di ata (z)=D A (z) V (z) 2 , 



(3) 



(4) 



actually, multiplied by rj~ 2 in the notation of Uzan et al. 
(2004). This quantity is reduced to the standard angular di- 
ameter distance only when the DD relation is strictly valid 
(rj = 1). In order to quantify the 77 parameter, Uzan et 
al. (2004) fixed Da(z) by using the cosmic concordance 
model (Spergel et al. 2003) while for D% ata (z) they con- 
sidered the 18 galaxy clusters from the Reese et al. (2002) 
sample for which a spherically symmetric cluster geometry 
has been assumed. By assuming 7/ constant, their statistical 
analysis provided rj = 0.91±o o 4 (1<t), and is therefore only 
marginally consistent with the standard result. 

Basset & Kunz (2004) used supernovae la (SNe la) data 
as measurements of the luminosity distance and the esti- 
mated Da from FRIIb radio galaxies (Daly & Djorgovski 
2003) and ultra compact radio sources (Gurvitz 1994, 1999; 
Lima & Alcaniz 2000, 2002; Santos & Lima 2008) in order 
to test possible new physics based on the following gener- 
alization of Eq. (JT|) 



Dl(z) 



D A {z){l + zf 



(1+*) 



exp 



dz' 



(5) 



where E(z) = H(z)/Hq is the dimensionless Hubble expan- 
sion, a quantity normalized to unity today. Note that for 
arbitrary values of a, the strict validity of the DD relation 
corresponds to (/3, 7) = (1, 0). By marginalizing on Qm, 
and Hubble parameters, Basset & Kunz (2004) found a 2c 
violation caused by excess brightening of SNIa at z > 0.5, 
perhaps owing to a lensing magnification bias. 

On the other hand, De Bernardis, Giusarma & 
Melchiorri (2006) also confronted the angular distances 
from galaxy clusters with luminosity distance data from 
SNe la to obtain a model-independent test. In order to 
compare the data sets they considered the weighted average 
of the data in seven bins and found that 77 = 1 is consis- 
tent on a 68% confidence level (la). However, one needs 
to be careful when using the Sunyaev-Zel'dovich effect to- 
gether with X-ray techniques for measuring angular diame- 
ter distances to test the DD relation because this technique 
also depends on its validity. Indeed, when the relation does 
not hold, Da(z) determined from observations is in general 
Dj^ ta (z) — DA(z)n 2 , which is reduced to Da only if r\ = 1. 
This means that De Bernardis and co-workers did not re- 
ally test the DD relation, at least not in a consistent way. 
In addition, these authors binned their data, and, as such, 
their results may have been influenced by the particular 
choice of redshift binning. 

Avgoustidis et al. (2010) also adopted an extended DD 
relation, D L = D A (l + z) 2+e , in the context of a flat ACDM 
model for constraining the cosmic opacity. The recent SN 
type la data compilation (Kowalski et al. 2008) was com- 
bined with the latest measurements of the Hubble expan- 
sion at redshifts in the range < z < 2 (Stern et al. 2010). 
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They found e 



-0.04] 



0.07 



(2a). It should be stressed, 



however, that the main goal of the quoted studies was to 
merely test the consistency between the assumed cosmo- 
logical model and the results provided by a chosen set of 
astrophysical phenomena. 

Below we explore a different route to test the DD re- 
lation by using two complementary, but independent, ap- 
proaches. First, we take the DD's validity for granted in 
order to access the galaxy cluster morphology. The ba- 
sic idea is very simple and may be described as follows. 
The usually assumed spherical geometry of clusters has 
been severely questioned after some analyses based on data 
from the XMM-Newton and Chandra satellites, which sug- 
gested that clusters are supposed to exhibit an elliptical 
surface brightness. In this context, by assuming the ACDM 
framework (WMAP7), we discuss the constraints coming 
from the validity of DD relation on the local physics, that 
is, when different assumptions about the cluster geometry 
are considered. In order to do that, we considered three 
samples of angular diameter distances from galaxy clus- 
ters obtained through Sunyaev-Zel'dovich effect and X- 
ray measurements. These samples differ by the assump- 
tions adopted to describe the clusters: (i) isothermal el- 
liptical, (ii) isothermal spherical (3 models (Dc Filippis et 
al. 2005), and non-isothermal spherical double (3 model 
(Bonamcntc et al. 2006). Secondly, we propose a consis- 
tent cosmological-model-independent test for the DD rela- 
tion by using subsamples of SNe la carefully chosen from 
Constitution data (Hicken et al. 2009) and the angular di- 
ameter distances from galaxy clusters. These topics were 
partially discussed by us (Holanda, Lima & Ribeiro 2010, 
2011) without considering the second possibility (isother- 
mal spherical (3 model) that has also been analyzed by De 
Filippis et al. (2005). Both approaches were separately in- 
vestigated, however, by avoiding all details of the Sunyaev- 
Zel'dovich effect and X-ray cluster physics. 

In this article, we intend to close the above described 
gaps by providing a closer discussion of the physics in- 
volved, and, for completeness, we also included the spheri- 
cal (3 model case. In addition, the influence of the different 
SNe la light-curve fitters on the model-independent test in- 
volving the Sunyaev-Zel'dovich effect, X-ray and SNe la is 
discussed. Finally, a joint analysis involving the different 
approaches is also investigated. The present study (based 
on complementary tests) suggests that clusters are endowed 
with an elliptical geometry as assumed by De Filippis et al. 
(2005), once the strict validity of DD relation is taken for 
granted. 

The paper is organized as follows. Sunyaev-Zel'dovich 
effect and X-ray surface brightness observations as a test 
for the DD relation are explored in <|2] The galaxy cluster 
samples used in this paper are presented in <J3J The con- 
sistence between the validity of the DD relation and the 
different assumptions used to describe the galaxy clusters 
usually adopted in the literature are discussed in fJU In 
Sj5] we discuss a new and model-independent cosmological 
test for the DD relation involving luminosity distances from 
SNe la and Da(z) from galaxy clusters. In §6 we study the 
influence of the different SNe la light-curve fitters on the 
results of the previous section, while the joint analysis is 
presented in §7. Finally, the main conclusions and future 
prospects are summarized in ^H] 



2. SZE/X-ray technique and the distance-duality 
relation 

An important phenomenon occurring in galaxy clusters is 
the Sunyaev-Zel'dovich effect (SZE), a small distortion of 
the CMBR spectrum provoked by the inverse Compton 
scattering of the CMBR photons passing through a pop- 
ulation of hot electrons. The SZE is proportional to the 
electron pressure integrated along the l.o.s., i.e., to the first 
power of the plasma density. The measured temperature 
decrement ATsz of the CMBR is given by (De Filippis et 
al. 2005) 



■ = f{v,T e ) 



p. To dl , 



(0) 



l.o.s. 



where T e is the temperature of the intra-cluster medium, 
k-Q the Boltzmann constant, T C mbr = 2.728°K is the tem- 
perature of the CMBR, or the Thompson cross section, 
m e the electron mass and f(v,T c ) accounts for frequency 
shift and relativistic corrections (Itoh, Kohyama & Nozawa 
1998; Nozawa, Itoh & Kohyama 1998). 

Other important physical phenomena occurring in the 
intra-galaxy cluster medium are the X-ray emission caused 
by thermal bremsstrahlung and line radiation resulting 
from electron-ion collisions. The X-ray surface brightness 
Sx is proportional to the integral along the line of sight 
of the square of the electron density. This quantity may be 
written as follows 



Sx — 



n 2 

A 

4tt£>£ 



n e A eH dl, 



(7) 



where A e # is the X-ray cooling function of the intra-cluster 
medium (measured in the cluster rest frame) and n e is the 
electron number density. It thus follows that the SZE and 
X-ray emission both depend on the properties (n e , T e ) of 
the intra cluster medium. 

As is well known, it is possible to obtain the angular 
diameter distance from galaxy clusters by their SZE and X- 
ray surface brightness observations. The calculation begins 
by constructing a model for the cluster gas distribution. 
Assuming, for instance, the spherical isothermal /3-model 
such that n e is given by (Cavaliere & Fusco-Fermiano 1978) 



n ^ r ) = ( 1 + ^ 



-3/3/2 



(8) 



equations (6) and (7) can be integrated. Here r c is the core 
radius of the galaxy cluster. This (3 model is based on the 
hydrostatic equilibrium equation and constant temperature 
(Sarazin 1988). In this way, we may write for the SZE 



AT SZ = AT [1 + ^7 



2X 1/2-3/3/2 



(9) 



where 9 C — r c /DA is the angular core radius and ATo is 
the central temperature decrement that includes all physi- 
cal constants and terms resulting from the line-of-sight in- 
tegration. More precisely: 

J rri 

AT = T CMBR f{v,T c f T B 2 c n e ^9 c D A g {(3/2) , (10) 



with 

. . r [3a- 1/2] 
9(a) = 



T 3a 



(11) 
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where T(a) is the gamma function and the others constants 
are the usual physical quantities. For X-ray surface bright- 
ness, we have 



S 



x 



~>X0 



1/2-3/3 



(12) 



where the central surface brightness Sxo reads 

D\A e H Pe/p-H 2 



S 



xo 



D|40F 



-ni 6 c D A g(fi). 



Here p is the molecular weight given by pi = p/riim p . 

One can solve equations (10) and (13) for the angular di- 
ameter distance by eliminating n e o and taking for granted 
the validity of DD relation. However, a more general re- 
sult appears when the DD relation is not regarded as being 
strictly valid. In this case one obtains 



D A = 



ATp 2 /m c c 2 y g(0) 
Sxo \k B T e oJ <?(/3/2) 2 #c 



4^/2/(^,^)2 (T< 



■ !Br) 2 Ct 2 



(1 + z c ) 4 r?(z) 2 



t~\ data _ — 2 
U A V ) 



(14) 



where z c is the galaxy cluster redshift. Therefore, as previ- 
ously stressed by Uzan et al. (2004), galaxy cluster observa- 
tions do not provide the angular diameter distance directly. 
In principle, instead of the real angular diameter distance, 
the measured quantity is D A lata {z) = Da(z) rj 2 (z). 

To proceed with our analysis, the quantity r)(z) defining 
the deformed DD relation (see Eq. ©) is parametrically 
described by two different expressions (Holanda, Lima & 
Ribeiro 2010; 2011) 



i) i](z) = 1 + rj z, ii) t](z) = 1 + r]oz/(l + z). 



(15) 



The first expression is a continuous and smooth one- 
parameter linear expansion, whereas the second one in- 
cludes a possible epoch- dependent correction that avoids 
the divergence at extremely high z. Naturally, one may ar- 
gue that these relations were not derived from first princi- 
ples. However, we stress that these expressions are very sim- 
ple and have several advantages such as a manageable one- 
dimensional phase space and a good sensitivity to obser- 
vational data. Clearly, the second parametrization can also 
be rewritten as rj(z) = 1 + 770(1 — a), where a(z) — (l + z) -1 
is the cosmic-scale factor. It represents an improvement 
with respect to the linear parametrization, since the DD 
relation becomes bounded regardless of the redshift values. 
Potentially, it will become more useful once higher redshift 
cluster data are made available. 

The above parametrizations are clearly inspired by sim- 
ilar expressions for the w(z)-equation of state parame- 
ter in time-varying dark energy models (see, for instance, 
Padmanabhan & Choudury 2003; Linder 2003; Cunha, 
Marassi & Santos 2007; Silva, Alcaniz & Lima 2007). In 
the limit of very low redshifts (z << 1), we have 77 = 1 and 
Di = Da as should be expected, and, more important for 
our subsequent analysis, the value 770 = must be favored 
by the Etherington principle. In other words, for a given 
data set, the likelihood of 770 must be peaked at 770 = 0, 
in order to satisfy the Etherington theorem. It should be 



remarked that Gongalves, Holanda & Alcaniz (2011) also 
adopted these expressions to explored the DD relation by 
using observations of gas mass fractions of galaxy clusters, 
whereas Lima, Cunha & Zanchin (2011) showed how to de- 
rive them from first principles based on possible theoretical 
modifications of the luminosity distance (without refraction 
effects). 



(13) 3. Galaxy cluster samples 



Below, the physical constraints encoded in the possibility 
of a deformed DD relation is explored by considering three 
samples of angular diameter distances from galaxy clus- 
ters obtained by combining their SZE and X-ray surface 
brightness observations. The first is defined by 38 angu- 
lar diameter distances from galaxy clusters in the redshift 
range 0.14 < z < 0.89, as given in the Bonamente et al. 
(2006) sample, where the hydrostatic equilibrium model 
and spherical symmetry was considered to describe the 
galaxy clusters. In order to construct a realistic model for 
the cluster gas distribution and include the possible pres- 
ence of the cooling flow, these authors modeled the gas den- 
sity with a function whose form is given below (Bonamente 
et al. 2006; La Roque et al. 2006), 



n e (r) = n eQ 



/ I 



9 X 



+ (1-/) 1 + 



3 3 



(16) 



This double /3-model for the density generalizes the sin- 
gle /3-model profile, introduced by Cavaliere and Fusco- 
Fermiano (1976) and the double f3 model proposed by Mohr 
et al. (1999). It has the freedom of following both the cen- 
tral spike in density and the more gentle outer distribution. 
The quantity n e o is the central density, / governs the frac- 
tional contributions of the narrow and broad components 
(0 < / < 1), r c i and r c i are the two core radii that describe 
the shape of the inner and outer portions of the density 
distribution and /3 determines the slope at large radii (the 
same (3 is used for both the central and outer distributions 
in order to reduce the total number of degrees of freedom) . 

On the other hand, the hydrostatic equilibrium and 
spherical symmetry hypotheses result in the condition 



dP 
dr 



(17) 



where P is the gas pressure, p g is the gas density and <p = 
—GM(r)/r is the gravitational potential due to dark matter 
and the plasma. Using the ideal gas equation of state for 
the diffuse intracluster plasma, P — p g ksT ' j \i m p , where p 
is the mean molecular weight and m p is the proton mass, 
one obtains a relationship between the cluster temperature 
and the cluster mass distribution, 



dT 
dr 



( pm p d(j) T dp g 
\ ks dr p g dr 
f pm p GM ^ T dp g 
V k B r 2 p„ dr 



(18) 
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Bonamente et al. (2006) combined hydrostatic equilibrium 
equations with a dark matter density distribution from 
Navarro, Frenk and White (1997), 



pDu{r) =N 



1 



(r/r s ){l + r/r s ) 2 



(19) 



where Af is a density normalization constant and r s is 
a scale radius. The parameters of these equations (n e o, 
/j fci, Tc2, /3, A/" and r s ) were obtained by calculating 
the joint likelihood C of the X-ray and SZE data in a 
Markov chain Monte Carlo method (Bonamente et al. 
2004). Summarizing, the cluster plasma and dark matter 
distributions were analyzed assuming hydrostatic equilib- 
rium model and spherical symmetry, thereby accounting for 
radial variations in density, temperature and abundance. 

The second sample is formed by 25 galaxy clusters in 
the redshift range 0.023 < z < 0.8 compiled by De Filippis 
et al. (2005). These authors re-analyzed archival X-ray data 
of the XMM-Newton and Chandra satellites of two samples 
for which combined X-ray and SZE analysis have already 
been reported using an isothermal spherical /3-model. One 
of the samples, compiled previously by Reese et al. (2002), is 
a selection of 18 galaxy clusters distributed over the redshift 
interval 0.14 < z < 0.8. The other one, the sample of Mason 
et al. (2001), has seven clusters from the X-ray limited flux 
sample of Ebeling et al. (1996). In this way, De Filippis et al. 
(2005) used an isothermal elliptical /3-model and an isother- 
mal spherical /3 model to obtain Da(z) measurements for 
these galaxy clusters samples. As discussed by De Filippis 
et al. (2005), the choice of circular rather than elliptical /3 
model does not affect the resulting central surface bright- 
ness or Sunyaev-Zeldovich decrement, the slope /3 differs 
slightly between these models, but significantly different 
values for core radius are obtained. The result was that the 
core radius of the elliptical /3-model is bigger than that of 
the spherical /3 model (see Fig. 1 in their paper). In a first 
approximation it was found that 6 e u = —£Z2i 9 C i rc , where 

i-~\~Gproj 

e pro j is the axial ratio of the major to the minor axes of the 
projected isophotes. Since Da oc angular diameter dis- 
tances obtained by using an isothermal spherical /3-model 
are overestimated compared with those from the elliptical 
/3-model. 

For the isothermal elliptical /3-model De Filippis et al. 
(2005) used a general triaxial morphology to describe the 
intra-cluster medium. They obtained 



sz 



where 



AT 



+ (eproj) V 



y c,proj 



1/2-3)9/2 



(20) 



AT = T CMBR f(u, T c ) °" T fcB J e n t -, , \ 



m c c 
Da 0c, P roj fgTeg 



/i 3 / 4 



(21) 



'Proj 



g(a) is given by equation (ITT1) . Da is the angular diameter 
distance to the cluster, 0i = x ii0 b s /£>A is the projected 
angular position (on the plane of the sky) of the intrinsic 
orthogonal coordinate Xi. h s , h is a function of the cluster 
shape and orientation, e pro j is the axial ratio of the major 
to the minor axes of the observed projected isophotes and 



#c,proj is the projection on the plane of the sky (p.o.s.) of 
the intrinsic angular core radius. 

Similarly, the X-ray surface brightness Sxo can be writ- 
ten as follows, 



~>x = Sxo 



2 Q 2 



+ ( e proj) ®2 



1/2-3/3 



($c,proj) 

where the central surface brightness Sxo reads 



(22) 



c _ ^eHD 2 A (pe/ Hh) 2 -DA ^c,proj 



eie 2 



9(P), (23) 



and fi is the molecular weight, given by fii = p/nim p . 

By eliminating n G o in the equations above and assum- 
ing the DD relation as valid, i.e., D A iata (z) = Da{z), De 
Filippis et al. (2005) obtained the observational quantity 
as written below 



D A (z) 



9{P) 



ATp 2 

Sxo \k B T eQ J ff (/3/ 2 ) : ^ proj 

A eg Q (Me/Mgj 

4^/2 /(l/;Te )2 ( TcMBR ) 2 a T 2 (l + z c y 



(24) 



The slope /3 of the profile and the projected core radius 
6c,proj were obtained by fitting the cluster surface bright- 
ness with an elliptical 2-D j3 model. For the isothermal 
spherical /3 model description, De Filippis et al. (2005) con- 
sidered the usual Eq. (8) and obtained angular distances 
with Eq. (14) and r] = 1. 

In Fig. 1 we plot the galaxy cluster samples. In Fig. la 
the filled circles (blue) and open squares (black) with the 
associated error bars (only statistical errors) stand for the 
De Filippis et al. (2005) isothermal elliptical /3 model and 
isothermal spherical /3 model, respectively. In Fig. lb we 
show the sample of Bonamente et al. (2006), where a non- 
isothermal spherical double (3 model was used to describe 
the galaxy clusters (red filled squares). 



4. Obtaining the shape of galaxy clusters by using 
the DD relation as constraint 

Many studies about the intra-cluster gas and dark mat- 
ter (DM) distribution in galaxy clusters have been limited 
to the standard spherical geometry (Reiprich & Boringer 
2002; Bonamente et al. 2006; Shang, Haiman & Verdi 2009). 
However, in the past few years observations of galaxy clus- 
ters based on Chandra and XMM data have shown that 
generally clusters exhibit elliptical surface brightness maps. 
Simulations have also predicted that DM halos show axis 
ratios typically on the order of w 0.8 (Wang & White 
2009), thereby disproving the spherical geometry assump- 
tion. In this line, the first determination of the intrinsic 
three-dimensional (3D) shapes of galaxy clusters was pre- 
sented by Morandi, Pedersen & Limousin (2010) by com- 
bining X-ray, weak-lensing and strong-lensing observations. 
Their methodology was applied to the galaxy cluster MACS 
J1423. 8+2404 and they found a triaxial galaxy cluster ge- 
ometry with DM halo axial ratios 1.53±0.15 and 1.44±0.07 
on the plane of the sky and along the line of sight, respec- 
tively. 
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Bearing in mind these results, we propose a new method 
to access the galaxy cluster morphology by taking the va- 
lidity of the DD relation as a constraint. The idea is very 
simple. Beacuse the samples shown in SJ3] were compiled 
under different geometric assumptions, we confront these 
underlying hypotheses with the validity of the DD relation. 
In principle, this kind of result provides an interesting ex- 
ample of how a cosmological (global) condition correlates 
to the local physics. In the application, one should also 
keep in mind that a deformed DD relation as given by Eq. 
3 naturally induces a more general result for the angular 
diameter distance from galaxy clusters via SZE and X-ray 
technique (see Eq. 14) 



D aata {z) = Da{z) ^ 



(25) 



In this line, our aim is to estimate the r/o parameter for 
each galaxy cluster sample for both parameterizations of 
■q(z) as given by Eq. (fT5")) . In our analyses (Holanda, Lima 
& Ribeiro 2011), Da(z) for each galaxy clusters is ob- 
tained from the WMAP (seven years) results by fixing 
the conventional flat ACDM model. The parameters of the 
simplest six-parameter ACDM model were recently deter- 
mined by Komatsu et al. (2011) by using the combination 
of seven-year data from WMAP with the latest distance 
measurements from baryon acoustic oscillations (BAO) in 
the distribution of galaxies (Percival et al. 2010) and the 
Hubble parameter Hq measurement presented by Riess et 
al. (2009). The basic results are n A = 0.725 ± 0.016 and 
h = 0.702 ± 0.014, and no convincing evidence for devia- 
tions from the minimal cosmic concordance model has been 
established. 

The theoretical angular diameter distance can be writ- 
ten as (Lima et al. 2003; Cunha et al. 2007) 



„ , , , 3000/1" 1 f* dz' 



(26) 



where h = Hq/1Q{) km s 1 Mpc 1 and the dimensionless 
function H(z'; f2 m ) is given by 



u= [n m (i + z') 3 + (i-n m )] 



1/2 



(27) 



Now, in order to constrain 770, let us evaluate the likeli- 
hood distribution function, e _x / 2 , where 



, ^{[r)(z)} 2 -[Vo b s(z)?Y 

x 2 = E — 



cr 



Vobs 



with [rf obs (z)] = D A lata {z)/D A {z) and 



'fobs 



D A {z) 



data 



T^data 
U A 



(28) 



(29) 



where o% MAP is the error in Da{z) associated to cosmo- 
logical parameters. The common statistical error contribu- 
tions for galaxy clusters are SZE point sources ±8%, X-ray 
background ±2%, Galactic N H < ±1%, ±15% for clus- 
ter asphericity, ±8% kinetic SZ and for CMBR anisotropy 
< ±2%. On the other hand, the estimates for system- 
atic effects are SZ calibration ±8%, X-ray flux calibra- 
tion ±5%, radio halos +3%, and X-ray temperature cal- 
ibration ±7.5%. Indeed, one may show that typical sta- 
tistical errors amount to nearly 20%, in agreement with 



other works (Mason et al. 2001; Reese et al. 2002, 2004), 
while for systematics we also find typical errors around + 
12.4% and - 12% (see also Table 3 in Bonamente et al. 
2006). In the present analysis we have combined the sta- 
tistical and systematic errors in quadrature for the galaxy 



2 

data 



°ltat + o-Lrf)- We note that in our X 



clusters (cr 

statistical analysis the asymmetric uncertainties present in 
the Bonamente et al. (2006) and De Filippis et al. (2005) 
samples were symmetrized by using the D'Agostini (2004) 
method. 

In Figs. 2a and 2b we plot the likelihood distribution 
function for the De Filippis et al. (2005) and Bonamente et 
al. (2006) samples. The 770 values and their errors (statis- 
tical + systematic errors) of our analysis are given below 
(see Table [J). 

For an isothermal elliptical (3 model, we see that the an- 
gular diameter distances from the De Filippis et al. (2005) 
sample provide an excellent fit and agree with the DD rela- 
tion at lcr confidence level. On the other hand, although 
it agrees with the DD relation at 2cr c.L, the analysis 
based on the isothermal spherical j3 model leads to 770 val- 
ues higher than those from the elliptical (3 case. Since in 
this case 770 > 0, the departures of the DD relation va- 
lidity indicates that the estimated angular distances with 
the spherical f3 model are overestimated with respect to 
those from the conventional flat ACDM (WMAP7). For 
the Bonamente et al. (2006) sample (non-isothermal spher- 
ical double f3 model) we also see that the strict validity of 
the reciprocity relation is only marginally compatible. The 
relative situation is not modified even when only clusters 
with z > 0.1 are considered in the De Filippis et al. (2005) 
sample. In such a situation, we obtain 770 = — 0.0441° 1 
(x 2 /d.o.f. — 15.9/17) for the linear parametrization, and 
?7o = -0.07ig;^ (x 2 /d-o.f. = 15.6/17) within la in the 



+0.11 
-0.1 



non-linear case for elliptical description and 770 = 0.186 
(x 2 /d.o.f. = 10.9/17) for the linear parametrization, and 
770 = 0.274t°;^ ( x 2 /d.o.f. = 10.7/17) within lcr in the 
non-linear case for spherical description. 

Therefore, we found no evidence for a distance-duality 
violation for the elliptical De Filippis et al. sample (2005). 
However, the same kind of analysis contradicts the spheri- 
cal symmetry hypothesis assumed in the Bonamente et al. 
(2006) sample and in the De Filippis et al. (2005) sample 
when a spherical geometry is assumed. These results are 
very interesting since they show how important the choice 
of geometry is in describing the clusters to obtain their dis- 
tances through SZE + X-ray measurements. We also see 
that the non-isothermal assumption of Bonamente et al. 
(2006) in their spherical description was not sufficient to 
satisfy the validity of the DD relation in the ACDM frame- 
work (WMAP 7, Komatsu et al. 2011). 

We recall that Bonamente et al. (2006) used three differ- 
ent models to describe the same 38 galaxy clusters: (i) the 
non-isothermal spherical double /3 model already discussed, 
(ii) an isothermal spherical (3 model, and (iii) the isothermal 
spherical (3 model excluding the central 100 kpc from the 
X-ray data. For the sake of completeness, we also obtained 
the ?7o values for the last two descriptions. For an isothermal 
spherical j3 model with linear and non-linear parametriza- 
tions we found r/ = — 0.09±g;g| ( X 2 /d.o.f. = 48.48/37) 
and 770 = — 0.13±g;gl ( X 2 /d.o.f. = 47.56/37), respectively. 
In the isothermal spherical f3 model excluding the central 
100 kpc from the X-ray data we obtained 770 = —0.14^5 07 
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Table 1. The 770 values and their errors (statistical + systematic errors) of our first analysis (section|?]). 



Isothermal elliptical /3 model x ACDM (WMAP7) 




770 (Linear case) = —0.056 ± —0.10 (la) 
7/0 (Non-linear case) = -0.088 ±0.14 (la) 


23.52/24 
22.56/24 


Isothermal spherical /3 model x ACDM (WMAP7) 


X'fd.o.f 


770 (Linear case) = 0.19 ± 0.12 (la) 
770 (Non- linear case) = 0.28 ± 0.18 (Iff) 


20.16/24 
19.44/24 


Non-isothermal spherical f3 model x ACDM (WMAP7) 


x'/d-o.f 


770 (Linear case) = —0.11 ± 0.06 (Iff) 
770 (Non-linear case) = -0.16 ± 0.08 (la) 


29.6/37 
29.23/37 



(x 2 /d.o.f. — 47.81/37) for the linear parametrization, and 
770 = — 0.2±g:?? (x 2 /d.o.f. = 48.40/37) in the non-linear 
case. All these determinations are at a la confidence level. 

Furthermore, Bonamente et al. (2006) also determined 
the Hubble constant Hq for each model. By assuming a flat 
ACDM model (CIm = 0.3 and Qa=0.7), the correspond- 
ing values at la c.l. are Hq = 76. 9+ 3 ' 4 , Hq 



and H = 77.6_ 4 ; 3 . More recently, Cunha, Marassi & Lima 
(2007) obtained Hq — 74±3.5 through a joint analysis with 
BAO by using the De Filippis et al. sample (2005), thereby 
alleviating the tension between between SZ + X-ray tech- 
nique and the CMB + BAO determination of Hq. However, 
since the Hubble constants obtained from all these models 
agree at la, one may conclude that the Hq value is not use- 
ful to decide which galaxy cluster model is more realistic. 



5. Testing the DD relation with galaxy clusters and 
SNe la 

In this section we apply a model-independent cosmological 
test for the DD relation based on the three samples de- 
fined in section 3. For Dl we considered two subsamplcs 
of SNe type la taken from Tables I and II of the Hicken et 
al. (2009) (Constitution data ), whereas values for Da are 
provided by the three samples of galaxy clusters discussed 
above. The SNe la redshifts of each subsample were care- 
fully chosen to coincide with those of the associated galaxy 
cluster sample (Az < 0.005), thereby allowing a direct test 
of the DD relation. For a given pair of data sets (SNe la, 
galaxy clusters), one should expect a likelihood of 770 peaked 
at 770 = 0, in order to satisfy the DD relation. Moreover, 
in our approach the data do not need to be binned as as- 
sumed in some analyses involving the DD relation (see, for 
instance, De Bernardis, Giusarma & Melchiorri 2006). 

from 



In Fig. 3a we plot Da multiplied by (1 + z ( 



luste 



the galaxy clusters sample compiled from the De Filippis 
et al. sample (2005) and Dl from our first SNe la subsam- 
ple. In Fig. 3b we plot the subtraction of redshift between 
clusters and SNe la. We see that the biggest difference is 
Az S3 0.01 for three clusters (open squares), while for the 
remaining 22 clusters we have Az < 0.005. In order to avoid 
the corresponding bias, we removed the three clusters from 
all subsequent analyses so that Az < 0.005 for all pairs. 

Similarly, in Fig. 4a we plot Da multiplied by (1 + 
Zduster) 2 , but now for the Bonamente et al. sample (2006) 
and Dl from our second SNe la subsample. In Fig. 4b we 
display the redshift subtraction between clusters and SNe 
la. Again, we see that for 35 clusters Az < 0.005. The 



biggest difference is again Az « 0.01 for three clusters, and, 
for consistency, these were also removed from our analysis. 

Let us now estimate the parameter 770 for each sam- 
ple and the two parametrizations defined by Eq. Q. We 
recall that in general the SZE + X-ray surface bright- 
ness observations technique does not produce Da(z), but 
D A iata (z) — Da(z) rj 2 . Consequently, if one wishes to test 
Eq. (1) with SZE + X-ray observations from galaxy clus- 
ters, the angular diameter distance must be replaced by 
D A iata (z) rf 2 in Eq. (3). In this way, we obtained r](z) = 
Di ata (z)(l + z cluster ) 2 /D L (z). 

Following standard procedure, the likelihood estimator 
is determined by a \ 2 statistics 



E 



[v(z) - Vobs(z)Y 



(30) 



Vobs 



where r] obs (z) = (1 + z cluster ) 2 D ^ ata (z) / 'D L (z) and 



(1 + Z c l us ter) 2 D d A ata 
(1 H~ ^cluster) 



'D, 



Di 



a nda 



(31) 



As in the previous section, we combined the statistical and 
systematic errors in quadrature for the angular diameter 
distance from galaxy clusters, and, as we remarked above, 
the asymmetric error bars were treated by the D'Agostini 
(2004) method. 

On the other hand, after nearly 500 discovered SNe la, 
the constraints on the cosmic parameters from luminosity 
distance are now limited by systematics rather than by sta- 
tistical errors. In principle, there are two main sources of 
systematic uncertainties in SNe la observations, which are 
closely related to photometry and possible corrections for 
light-curve shape (Hicken at al. 2009). However, at the mo- 
ment it is not so clear how one is to estimate the overall 
systematic effects for these standard candles (Komatsu et 
al. 2011), and, therefore, we exclude them from the follow- 
ing analysis. The basic reason is that systematic effects from 
galaxy clusters seem to be larger than those of SNe observa- 
tions, but their inclusion does not affect the results validity 
of the distance-duality relation very much. In this way, fol- 
lowing Holanda, Lima & Ribeiro (2010), we included only 
statistical errors of the SNe la magnitude measurements. 

In Figs. 5a and 5b we plot the likelihood distribution 
function for each sample. The 770 values with their errors 
(statistical and systematic errors) are provided in Table 
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[21 Clearly, the confrontation between the angular diame- 
ter distances from the elliptical De Filippis et al. (2005) 
sample with SNc la data disagree moderately with the 
reciprocity relation (the DD relation is marginally satis- 
fied in 2a) and the 770 values are mostly negative. Since 
{l + z c iuster) 2 D^ ata {z)/D L {z) (x 1 + t] z, negative r] values 
indicate luminosity distances overestimated with respect to 
the angular diameter ones. This tension between Sne la and 
the elliptical De Filippis et al. (2005) sample arises because 
SNe la samples prefer universes with a higher J7a value 
than those of the WMAP7 results. Indirectly, our analysis 
support the tension found by Wei (2010), who used SNe la 
constitution data and observations of the CMBR anisotropy 
plus BAO separately to constrain the ui a parameter of the 
dark energy equation of state given by u> = cuo+ui a z/(l + z). 
The SNe la constitution data indicated a strong negative ui a 
(ui a — 11) and the phantom energy divide line (uj < — 1) 
could be crossed at recent redshifts. Furthermore, the con- 
stitution data disagreed not only with the CMB and BAO 
observations, but also with other SNe la datasets such as 
Davis (Davis et al. 2007) and SNLS (Astier et al. 2006). 
On the other hand, when an isothermal spherical /3 model 
is used, the DD relation is satisfied in la. As discussed ear- 
lier, this concordance occurs because an isothermal spher- 
ical /? model yields angular distances overestimated com- 
pared to the elliptical model and WMAP7 results. For the 
Bonamente et al. (2006) sample, where a non-isotermal 
spherical double (3 model was assumed to describe the clus- 
ters, we see that the DD relation is not obeyed even at 
3a. We note that our results do not change if one replaces 
^cluster by zsnc, where zsNe is the redshift of SNe la. This 
shows that a Az < 0.005 is sufficient to implement our 
direct test. 

In this context, we recall that Li, Wu and Yu (2011) 
rediscussed this independent cosmological test using the 
latest Union2 SNe la data and the angular diameter dis- 
tances from galaxy clusters. However, they removed six and 
twelve data points, respectively, from the De Filippis et al. 
(2005) and Bonamente et al. (2006) samples and obtained 
a more serious violation of the DD relation. These authors 
also reexamined the DD relation by postulating two more 
general parameterization forms, namely 77(2) = 770 + 7712 
and 77(2) = r/o + 7/12/(1 + 2), and they found that con- 
sistencies between the observations and the DD relation 
are markedly improved for both samples of galaxy clusters. 
Nair, Jhingan and Jain (2011) also discussed the strict va- 
lidity of the DD relation using the latest Union2 SNe la 
data and the angular diameter distances from galaxy clus- 
ters, FRIIb radio galaxies and mock data. They proposed 
six different (one and two indices) parametrizations (in- 
cluding, as particular cases, those adopted by Holanda et 
al. 2010) in an attempt to determine a possible redshift 
variation of the DD relation. As physically expected, their 
results depend both on the specific parametrization and the 
data samples considered. In particular, they conclude that 
one index parameterization, namely i~jv — 7 7s/(l + z ) and 
Vv = % exp{[z/(l + z)]/(l + 2)}, do not support the DD 
relation for the given data set. 

More recently, Meng, Zhang and Zhan (2011) also rein- 
vestigated the model-independent test by comparing two 
different methods and several parametrizations (one and 
two indices) for 77(2). Their basic conclusion agrees with our 
studies in the sense that the triaxial ellipsoidal model is sug- 
gested by the model-independent test at la while the spher- 



ical j3 model can only be accommodated at a 3a confidence 
level. As we shall see, this uncertainty can also be robustly 
resolved by considering a joint analysis involving the two 
treatments discussed here (see section 7). In principle, since 
all results somewhat depend on the assumed 77(2) form, it is 
important to understand the effect of different parametriza- 
tions. However, we recall that at extremely low redshifts 
the DD relation reduces to unity since Dl(z) = Da{z), 
and, therefore, the constant parameter appearing in the 
proposed 77(2) expressions should be fixed to unity. 

At this point it is interesting to know the influence of 
SNe la light-curve fitters on the model-independent test of 
the distance duality relation and some of its consequences 
for the local physics. This topic has attracted a lot of atten- 
tion in the recent literature (Kessler et al. 2009; Bengochea 
2010, Hao Wei 2011), and its connection with the DD rela- 
tion and galaxy cluster geometry it will be discussed next. 



6. SNe la light-curve fitters, DD relation and 
cluster geometries 

In the previous section, we tested the DD relation by using 
luminosity distances from SNe la based on the Constitution 
sample and the angular diameter distances from three 
galaxy clusters samples obtained by their SZE and X- 
rays surface brightness measurements. However, in the 
Constitution compilation different SNe la samples were an- 
alyzed by four light-curve fitters (SALT, SALT2, MLCS31 
and MLCS17) to test consistency and systematic differ- 
ences: 397 SNe la with SALT, 351 SNe la with SALT2, 
366 SNe la with MLCS31 and 372 SNe la with MLCS17. 
The four fitters were seen to be relatively consistent with 
the light-curve-shape and color parameters, but that still 
leaves room for improvement, and they provide a consider- 
able amount of systematic uncertainty to any analysis. In 
this way, it would be interesting to investigate the influence 
of different SNe la light-curve fitters on our previous results. 
Unfortunately, it was not possible to find the same SNe's 
la used here in section [5] in the Constitution sample, but 
analyzed with the four methods for a direct comparison. 

On the other hand, Kessler et al. (2009) presented the 
Hubble diagram for 103 SNe la with redshifts 0.04 < 2 < 
0.42, discovered during the first season of the Sloan Digital 
Sky Survey-II (SDSS-II). These data filled the redshift 
desert between low- and high-redshift SNe la of the pre- 
vious SNe la surveys. In addition, these authors included a 
comprehensive and consistent reanalysis of other datasets 
(ESSENCE, SNLS, HST), resulting in a combined sample 
of 288 SNe la. Reanalysis was performed by using two light- 
curve fitters: SALT2 (Guy et al. 2007) and MLCS2K2 (Jha 
et al. 2007). MLCS2K2 calibration uses a nearby train- 
ing set of SNe la assuming a close to linear Hubble law, 
while SALT2 uses the whole data set to calibrate empirical 
light curve parameters, and a cosmological model must be 
assumed in this method. Typically a ACDM or a wCDM 
model is assumed. Consequently, the SNe la distance mod- 
uli obtained with SALT2 fitter retain a degree of model de- 
pendence (Bengochea 2011). Kessler et al. (2009) combined 
these 288 SNe la with measurements of baryon acoustic os- 
cillations from the SDSS luminous red galaxy sample and 
with cosmic microwave background temperature anisotropy 
measurements from the WMAP for estimating the cos- 
mological parameters uj and SIm by assuming a spatially 
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Table 2. The (770) values and their errors (statistical + systematic errors) of our second analysis (section [5]) . 



Isothermal elliptical j3 model x SNe la 


X Z /d.o.f 


770 (Linear case) = -0.28 ± 0.21 (Iff) 
770 (Non-linear case) = -0.43 ± 0.29 (la) 


21.2/21 
21/21 


Isothermal spherical /3 model x SNe la 


X'fd.o.f 


770 (Linear case) = 0.14 ± 0.26 (Iff) 
770 (Non-linear case) = 0.20 ± 0.36 (Iff) 


17.01/21 
16.80/21 


Non- isothermal spherical /3 model x SNe la 


X'fd.o.f 


7/0 (Linear case) = —0.39 ±0.11 (Iff) 
770 (Non-linear case) = —0.61 ± 0.16 (la) 


30.26/34 
28.9/34 



flat cosmological model. The results from these two light- 
curve fitters disagreed: ui = — 0.76± 0.07(stat)±0.11(syst), 
n M = 0.307 ± 0.019(stat)±0.023(syst) using MLCS2K2 
and u) = -0.96 ± 0.06(stat)±0.12(syst), ft M = 0.265 ± 
0.016(stat)±0.025(syst) using SALT2. This discrepancy 
raised the question of which method is more reliable be- 
cause it is not possible to definitively determine from the 
current data that either method is better or incorrect. The 
overall conclusion was that the cosmological parameter ui 
lies between —1.1 and —0.7 (Kessler et al. 2009). Exploring 
nonstandard cosmology scenarios, Sollerman et al. (2009) 
found that more exotic models provided a better fit to the 
SNe la data than the ACDM model when the MLCS2K2 
light-curve fitter was used. However, when the SN la data 
were analyzed by using SALT2 light-curve fitter, the stan- 
dard cosmological constant model agreed better to those 
data. 

In this section, by using the SNe la data from SDSS-II 
(Kessler et al. 2009), we explore the influence of the two dif- 
ferent SNe la light-curve fitters on our model-independent 
cosmological test for the DD relation. For Dl we consid- 
ered two subsamples of SNe la taken from SDSS-II (2007) 
where, again, the SNe la redshifts of each subsample were 
carefully chosen to coincide with those of the associated 
galaxy cluster sample (Az < 0.005), allowing a direct test 
of the DD relation. Furthermore, in this case, each subsam- 
ple of SNe la was analyzed by the MLCS2K2 and SALT2 
light-curve fitters separately. For Da we used the angu- 
lar diameter distance samples of galaxy clusters discussed 
above: the isothermal elliptical and spherical /? model sam- 
ples from De Filippis et al. (2005) and the non-isothermal 
spherical double j3 model sample from Bonamente et al. 
(2006). 

In Fig. 6a we plot the subtraction of redshift between 
clusters and SNe la for the De Filippis et al. (2005) sam- 
ple. We see that the biggest difference is Az w 0.007 for 
three clusters (open squares), whereas for the remaining 22 
clusters we have Az < 0.005. In order to avoid the cor- 
responding bias, the three clusters were removed from all 
analyses presented here so that Az < 0.005 for all pairs. In 
Fig. 6b we display the redshift subtraction between clus- 
ters and SNe la for the Bonamente et al. (2006) sample. 
In this case we see that for 33 clusters Az < 0.005. The 
biggest difference is Az ps 0.015 for five clusters, and, for 
consistency, they were also removed from our analysis. 

In Figs. 7a and 7b we plot the likelihood distribution 
function for the elliptical and spherical De Filippis et al. 
(2005) samples, respectively. The results for Bonamente et 
al. (2006) sample are shown in Fig. 8. Following Kessler 



et al. (2009) and Sollerman et al. (2009), we added in our 
analyses an additional intrinsic dispersion of 0.16 mag to 
the uncertainties output by the MLCS2K2 light-curve fitter 
and 0.14 mag for the SALT2 light-curve fitter. In Table [3] 
we display our results, where the errors are in 2a (statistical 
plus systematic errors). 

Interestingly, we also obtained conflicting results be- 
tween the two light-curve fitters. In our analysis involv- 
ing the elliptical De Filippis et al. (2005) sample (for both 
light-curve fitters), the DD relation validity is obtained 
at la (c.l.) with the SALT2 method. However, it is only 
marginally compatible with MLC2K2 method because the 
770 values are considerably negative. This result points to 
an overestimated luminosity distances when SNe la are an- 
alyzed with the MLC2K2 method. This result is not sur- 
prising since the De Filippis et al. (2005) sample agrees 
with the DD relation validity in the context of the ACDM 
model (WMAP7) (see sectiong]) and the SALT2 light-curve 
fitter also favors the same cosmology according to the anal- 
yses of Kessler et al. (2009) and Sollerman et al. (2009). 
Curiously, the best verification of the DD relation valid- 
ity is obtained by confronting the spherical sample of De 
Filippis et al. (2005) with SNe la when the MLC2K2 light- 
curve fitter is adopted. This probably occurs because both 
distances are overestimated compared to distances from the 
elliptical model and the ACDM model (WMAP7). For the 
Bonamente et al. (2006) sample the DD relation validity 
is marginally compatible (2a) with SALT2 and at 3<r with 
MLCS2K2. 



7. Searching for the true geometry of clusters: a 
joint analysis 

In the previous sections we have discussed how the valid- 
ity of the DD relation can be used to fix the true cluster 
geometry by using two independent approaches: (i) The 
mo del- dependent test based on the SZE/X-ray technique 
in the context of the ACDM model (WMAP7 analysis), 
and (ii) The model-independent test using the combination 
of SZE/X-ray for galaxy clusters and SNe la data. For the 
former case the elliptical description provides the best fit 
supporting the standard DD relation in ACDM (Komatsu 
et al. 2011). On the other hand, the isothermal spherical ft 
model is more compatible with the validity of the DD re- 
lation in the model-independent test involving supernovas. 
Since likelihoods for both descriptions in both approaches 
(see figs 2a and 5a) are compatible with each other in at 
least 2-sigma, it is convenient to perform a joint analysis in- 
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Table 3. The 770 values for our third analysis (section [5]) . All errors are at la. 



Isothermal elliptical f3 model 


SALT2 




MLCS2K2 


X z /d-o.f 


770 (Linear case) 
770 (Non-linear case) 


-0.24 ±0.24 
-0.34 ± 0.34 


19.74/21 
19.53/21 


-0.43 ±0.21 
-0.63 ±0.32 


21.2/21 
20.8/21 


Isothermal spherical /3 model 


SALT2 




MLCS2K2 


x'/d-o.f 


770 (Linear case) 
770 (Non-linear case) 


0.27 ±0.27 
0.4 ±0.4 


16/21 
15.7/21 


0.01 ±0.34 
0.07 ±0.35 


15.7/21 
15.5/21 


Bonamente et al. Sample 


SALT2 


x'/d-o.f 


MLCS2K2 


X'ld.o.f 


770 (Linear case) 
770 (Non-linear case) 


-0.22 ±0.14 
-0.28 ±0.24 


21.76/32 
21.2/32 


-0.28 ±0.10 
-0.45 ±0.15 


24/32 
23.7/32 



volving these independent and complementary treatments 
in order to choose the better description. 

Naturally, one may argue that the joint analysis will be 
plagued by the underlying tension between SNe and CMB 
data. In this context, we recall that this tension was recently 
discussed by taking into account weak-lensing effects caused 
by inhomogcneitics at low and intermediate rcdshifts by 
Amendola et al. (2010. By using the Union data (2008), 
they showed that the inclusion of lensing moves the best-fit 
model significantly toward the flat LCDM model (owing to 
corrections induced by the lens convergence on the distance 
modulus). This treatment is beyond the scope of our paper. 
Nevertheless, this treatment is valid in the same sense as the 
early analysis involving SNe la and CMB were to the cosmic 
concordance model. As we shall see, the joint analysis below 
confirms that the elliptical isothermal model provides the 
best description for galaxy clusters. 

Now, by adding the \ 2 statistics for the different ap- 
proaches, the resulting likelihoods for the 770 parameter were 
again obtained. The main findings for both parametriza- 
tions (linear and non-linear cases) are summarized in Table 
4. Note that we restricted our attention to the De Filippis 
et al. (2005) samples (elliptical and spherical isothermal 
assumptions). 

From Table 4, we see that the elliptical geometry agrees 
better with the strict validity of the duality relation when 
both approaches are considered. This result remains valid 
regardless of the supernova sample (Constitution or SDSS) 
adopted, the parameterization (linear or non-linear) , or the 
kind of light fitter used. 

In Fig. 9, we plotted the likelihoods function for the 770 
parameter as obtained from our joint analysis by consider- 
ing again the De Filippis et al. (2005) samples (elliptical and 
spherical isothermal geometries). Fig. 9a is the joint analy- 
sis for the WMAP + SNe la by using the Constitution sam- 
ple while in Figs. 9b and 9c, we considered SDSS sample, 
but now taking into account the different light-curve fitters 
(SALT2 and MLCS2k2). The horizontal lines are cuts in the 
regions of 68.3 (Ict) per cent and 95.4 (2a) per cent proba- 
bility. Our joint analysis shows that an elliptical geometry 
is suggested by the existing data even when we take into 
account the SNe la light-curve fitters. Because several SZE 
surveys are in progress, our results call attention to a basic 
difficulty involving the spherical assumption to describe the 
cluster morphology when independent cosmological probes 
are considered. 



8. Conclusions 

We explored some physical consequences of a deformed 
distance-duality relation, 77(2) = Z?l(1 + z)^ 2 /Da, based 
on observations of galaxy clusters using SZE and X-ray 
surface brightness. The r](z) parameter was described by 
two distinct forms, 77 = 1 + 770-2 and 77 = 1 ± T}qz/(1 + z), 
thereby recovering the standard equality between distances 
for extremely low redshifts. 

Initially, we discussed the consistency between the strict 
validity of the distance-duality relation and the underlying 
assumptions about the galaxy cluster geometries. In our 
analysis we used angular diameter distances from galaxy 
clusters of the De Filippis et al. (2005) and Bonamente et 
al. (2006) samples. The former sample consists of 25 data 
sets where the galaxy clusters were described by an isother- 
mal elliptical j3 model and a spherical f3 model, whereas 
the latter sample consists of 38 data sets where a non- 
isothermal spherical double f3 model was used. We showed 
that the elliptical geometry is more consistent with no vi- 
olation of the distance duality relation in the context of 
ACDM (WMAP7 data). In the elliptical case of the De 
Filippis et al. (2005) sample (see Fig. 2a and table [1} we 
found 7/ = -0.0561"^ and 770 = -0.08812;" for linear and 
non-linear parametrizations at la (statistical plus system- 
atic errors), respectively. The analysis with the isothermal 
spherical /3 model led to higher and positive 770 values than 
those from elliptical j3 model, indicating that the estimated 
angular distances with the spherical f3 model are overesti- 
mated compared to those from ACDM (WMAP7). On the 
other hand, the non-isothermal spherical double ft model 
(see Fig. 2b) was only marginally compatible with 770 = 0, 
such as 770 = -Q.12±g;^ and 770 = — 0.175±g ;^| for lin- 
ear and non-linear parameterizations at 2a (statistical plus 
systematic errors) , respectively. In this way, our analysis re- 
vealed that the elliptical model is compatible with the du- 
ality relation at la, whereas the non-isothermal spherical 
model (in Bonamente et al. 2006 sample) is only marginally 
compatible at 3a in the ACDM framework. 

Furthermore, we discussed a new and model- 
independent cosmological test for the distance-duality 
relation. We considered the three angular diameter dis- 
tance samples from galaxy clusters, which are obtained 
by using SZE and X-ray surface brightness, together the 
luminosity distances given by two subsamples of SNe la 
taken from the Constitution data. For each subsample, 
the redshifts of the SNe la were carefully chosen to 
coincide with those of the associated galaxy cluster sample 
(Az < 0.005). Our results showed that the confrontation 
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Table 4. The 770 values for our joint analysis. All errors are at la and the x 2 /d.o.f. is given in parentheses. 



Isothermal elliptical j3 model 


WMAP+SNe la 
Constitution 


WMAP+SNe la 
SDSS(SALT2) 


WMAP+SNe la 
SDSS(MLCs2k2) 


770 (Linear case) 
770 (Non-linear case) 


-0.0961^(0.96) 
-0.14±°;^(0.97) 


-0.081-- (0.93) 
-0.121^(0.95) 


-O.lli;;;- (0.97) 
-0.177l°; 1 «(0.96) 


Isothermal spherical /3 model 


WMAP + SNe la 
Constitution 


WMAP+SNe la 
SDSS(SALT2) 


WMAP+SNe la 
SDSS(MLCS2k2) 


770 (Linear case) 
770 (Non-linear case) 


0.178^:^(0.76) 
0.263^:^(0.75) 


0.1951^(0.77) 
0.29lS;i*(0.76) 


0.16l u ;-(0.75) 
0.241^^(0.75) 



between the angular diameter distances from the De 
Filippis et al. (2005) sample (elliptical model) with SNe 
la data pointed to a moderate violation of the reciprocity 
relation (the DD relation was marginally satisfied in la) 
and the 770 values were predominantly negative. Since 
(1 + z) 2 D^ ata (z)/D L (z) oc 1 + 7702;, negative 770 values 
indicated luminosity distances overestimated with relation 
to angular diameter ones. This tension between Sne la and 
the elliptical De Filippis et al. (2005) sample arises because 
SNe la samples prefer universes with higher Q\ values 
than the WMAP 7 results. When an isothermal spherical 
/3 model is adopted, the DD relation had been satisfied 
at la with 770 values preferably positive. Seemingly, this 
concordance occurred only because an isothermal spherical 
P model yields angular distances overestimated in compar- 
ison to an elliptical model and WMAP7 results. For the 
Bonamente et al. (2006) sample, where a non-isotermal 
spherical double (3 model was assumed to describe the 
clusters, we saw that the DD relation is not obeyed even 
at 3a. For this case we obtained 770 = — 0.421^34 and 
770 = —0.66i ' 5 for linear and non-linear rj{z) param- 
eterizations in 3cr (statistical plus systematic errors), 
respectively. 

Moreover, by using the SNe la data from SDSS-II 
(Kessler et al. 2009), we explored the influence of the two 
different SNe la light-curve fitters in our model-independent 
cosmological test for the DD relation. In this way, we 
considered two subsamples of SNe la taken from SDSS-II 
(2007) for Dl where, again, the SNe la redshifts of each sub- 
sample were carefully chosen to coincide with those of the 
associated galaxy cluster sample (Az < 0.005). In this case, 
each subsample of SNe la was analyzed by the MLCS2K2 
and SALT2 light-curve fitters separately. Interestingly, we 
obtained conflicting results between the two light-curve fit- 
ters. For both light-curve fitter-method the best fits are ob- 
tain in an analysis involving the elliptical De Filippis et al. 
(2005) sample. In this case, the DD relation validity was ob- 
tained in la with the SALT2 method and it was marginally 
compatible with MLC2K2 method with 770 values consid- 
erably negative. This result points to overestimated lumi- 
nosity distances when SNe la are analyzed with MLC2K2 
method. For this light-curve fitter, the isothermal spherical 
(3 model provided that the DD relation was satisfied at la. 
For all tests of the Bonamente et al. (2006) sample, the DD 
relation validity is obtained marginally at 2a (SALT2) and 
at 3<T (MLCS2K2). We stress that our results are indepen- 
dent of the redshifts chosen to realize the test (of the galaxy 
clusters or SNe la), thereby showing the robustness of the 
method and supporting Az < 0.005 as a fair choice. 



Finally, we applied a joint analysis involving these in- 
dependent and complementary treatments by adding the 
X 2 statistics for the different approaches. Interestingly, we 
saw that the elliptical geometry agrees better with the 
strict validity of the duality relation when both approaches 
were considered. This result remains valid regardless of the 
adopted SNe la sample (Constitution or SDSS), the param- 
eterization (linear or non- linear) , or the kind of light-curve 
fitter used (see Table 4). The non-isothermal spherical dou- 
ble /3 model was only marginally compatible with the DD 
relation validity in all treatments. 

Summarizing, the statistical analysis presented here 
provides new evidence that the true geometry of clusters 
has an elliptical form. In principle, it is remarkable that a 
local property such as the geometry of galaxy clusters may 
be constrained by a global argument like the one provided 
by the cosmological distance-duality relation. This result 
also reinforces the interest in the observational search for 
SZE and X-ray from clusters at high redshifts. In the 
near future, when more and larger samples with smaller 
statistical and systematic uncertainties become available, 
the method proposed here, based on the validity of the 
distance-duality relation, can improve the limits on the 
possible cluster geometries and explore systematic errors 
in SNe la and galaxy clusters observations. 
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Fig. 1. a) Filled (blue) circles and open (black) squares with the associated error bars (only statistical errors) stand for 
the De Filippis et al. (2005) samples: isothermal elliptical (3 model and isothermal spherical j3 model, respectively, b) The 
sample of Bonamente et al. (2006) where a non-isothermal spherical (3 model was used to describe the galaxy clusters. 




Ho \ 



Fig. 2. a) Likelihood distribution functions for De Filippis et al. (2005) sample. The solid and the dotted blue lines are 
likelihood functions corresponding to linear and non-linear rj(z) parametrizations for the isothermal elliptical (3 model. The 
dashed and the dashed-dotted red lines are likelihood functions corresponding to linear and non-linear parametrizations 
for the isothermal spherical (3 model, b) The likelihood distribution functions for Bonamente et al. (2006) sample. The 
solid blue line and the dotted red line correspond to linear and non-linear parametrizations. 
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Fig. 3. a) Galaxy clusters and SNe la data. The open blue, open black, and filled red circles with the associated error bars 
stand for the galaxy clusters described with elliptical and spherical fi models of the De Filippis et al. (2005) (statistical 
+ systematical errors) and SNe la (only statistical errors) samples, respectively, b) The redshift subtraction for the same 
pair of cluster-SNe la samples. The open squares represent the pairs of points for which Az ps 0.01. 
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Fig. 4. a) Galaxy clusters and SNe la data. The open (blue) and filled (red) circles with the associated error bars stand for 
the Bonamente et al. (2006) (statistical + systematical errors) and SNe la (only statistical errors) samples, respectively, 
b) The redshift subtraction for the same pair of cluster-SNe la samples. As in Fig. 1 b, the open squares represent the 
pairs of points with the biggest difference in redshifts (Az ss 0.01). 




Fig. 5. a) Likelihood distribution functions for De Filippis et al. (2005) sample. The solid and the dotted blue lines are 
likelihood functions for linear and non-linear parametrizations corresponding to the isothermal elliptical /3 model. The 
dashed and the dashed-dotted red lines are likelihood functions for linear and non-linear parametrizations corresponding 
to the isothermal spherical (3 model, b) The likelihood distribution functions for Bonamente et al. (2006) sample. 
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Fig. 6. a) Redshift subtraction for the same pair of cluster-SNe la samples of De Filippis et al. (2005). The open squares 
represent the pairs of points for which Az > 0.005. b) The redshift subtraction for the same pair of cluster-SNe la samples 
for Bonamente et al. (2006). The open squares represent the pairs of points for which Az > 0.005. 
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Fig. 7. a) Likelihood distribution functions for the De Filippis et al. (2005) sample (elliptical (3 model) for both r)(z) 
parametrizations and both SNe la light-curve fitters, b) The likelihood distribution functions for the De Filippis et al. 
(2005) sample (spherical f3 model) for both r](z) parametrizations and both SNe la light-curve fitters. 
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Fig. 9. a) Joint analysis for the De Filippis et al. (2005) samples and WMAP + SNe la (Constitution Sample), b) Joint 
analysis for the De Filippis et al. (2005) samples and WMAP + SNe la (SDSS-SALT2). c) Joint analysis for the De 
Filippis et al. (2005) samples and WMAP + SNe la (SDSS-MLCS2k2). 
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